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Abstract: We present the two- loop QCD helicity amplitudes for quark-quark and 
Qh! quark-antiquark scattering. These amplitudes are relevant for next-to-next-to-leading or- 

p il der corrections to (polarized) jet production at hadron colliders. We give the results in the 

't Hooft-Veltman and four-dimensional helicity (FDH) variants of dimensional regulariza- 
tion and present the scheme dependence of the results. We verify that the finite remainder, 
after subtracting the divergences using Catani's formula, are in agreement with previous 
results. We also provide the amplitudes for gluino-gluino scattering in pure M = 1 super- 
5^ I symmetric Yang-Mills theory. We describe ambiguities in continuing the Dirac algebra to 

D dimensions, including ones which violate fermion helicity conservation. The finite re- 
mainders after subtracting the divergences using Catani's formula, which enter into physical 
quantities, are free of these ambiguities. We show that in the FDH scheme, for gluino-gluino 
scattering, the finite remainders satisfy the expected supersymmetry Ward identities. 
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1. Introduction 



In the coming years, our field is looking forward to the ongoing experiments at the Teva- 
tron at Fermilab and the future ones at the Large Hadron Collider at CERN for unlocking 
the physics of the electroweak symmetry breaking scale. At hadron colliders, for large 
momentum transfer, the most copious events are hadronic jets. To explore the validity of 
the Standard Model at the shortest possible distances, therefore, it would be helpful to 
determine jet production cross sections with high precision. Calculations of jet produc- 
tion at next-to-leading order (NLO) in the strong coupling constant as 0, agree well 
with the data over a broad range of transverse momentum. Still, the NLO predictions 
have an uncertainty from higher-order corrections, which is traditionally estimated using 
dependence on the renormalization and factorization scales, of order 10% or more. For 
very large momentum transfer the predictions can be improved by resumming threshold 
logarithms Q. There are also sizable uncertainties associated with the experimental in- 
put to the parton distribution functions Nevertheless, an exact next-to-next-to-leading 
order (NNLO) computation would be desirable. An important step has recently been ac- 
complished with the computation of the three-loop splitting function by Moch, Vermaseren 
and Vogt |p. There has also been some earlier work on global fits to the data Q within an 
approximate next-to-next-to-leading order (NNLO) framework |§. Once combined with 
the matrix elements in complete programs, this should considerably reduce the renormal- 
ization and factorization scale uncertainties in production rates. For a summary of the 
various expected improvements see, for example, ref. @. 

Recent years have seen rapid progress in our ability to compute two-loop matrix 
elements, especially when there is dependence on more than a single kinematic vari- 
able 0, |ll|, m, ll, m, m, m, ll, m, II], m, H. Much of tWs progress has rehed 
on new developments in loop integration |24, 25, 27, 28, 29, 30 1 and in understanding 



the infrared divergences of the theory |31]. 

An NNLO calculation of jet production requires requires six-point tree-level, one-loop 
five-point amplitudes and two-loop four-point amplitudes. The tree amplitudes for six 



external partons ||3^, ^] and the one-loop amplitudes for five external partons have 
been determined some time ago. Anastasiou, Glover, Oleari, and Tejeda-Yeomans have 
provided the NNLO interferences of the two-loop amplitudes with the tree amplitudes, for 
all QCD four-parton processes, summed over all external helicities and colors [15|. The 
helicity amplitudes for gg — > gg were presented in ref. ||l8|. The qq —>■ gg and qg — > 
qg helicity amplitudes were presented in refs. |21, 22 1. Recently, while preparing this 
paper. Glover presented the four-quark helicity amplitudes [^. Here we present the same 
amplitudes using somewhat different methods, as well as the M = 1 super symmetric version 
of these amplitudes. We also describe ambiguities in the amplitudes that first arise in four- 
fermion amplitudes at NNLO. 

For jet production in collisions of unpolarized hadrons, which is the main phenomeno- 
logical application of the amplitudes, the additional helicity and color information contained 
in the helicity amplitudes is not necessary. However, for the case of polarized proton scat- 
tering at the relativistic heavy ion collider (RHIC) at Brookhaven, the helicity amplitudes 
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are of direct relevance for improving predictions to NNLO accuracy. This may help with the 
determination of the poorly-known polarized gluon distribution in the proton ||35[1 , which 
is currently available through NLO p6[ |. 

Many formal properties of scattering amplitudes are simpler in a helicity basis. As a 
striking recent example, Witten has linked tree-level helicity amplitudes to a twistor space 
topological string theory |^] leading to simple [|^] and efficient |Q rules for dealing with 
general pO] tree amplitudes in massless gauge theories. It may also lead to new insights 



into loop calculations |41]. Other examples are supersymmetry Ward identities [^, 43 1, 
the behavior of amplitudes as momenta become coUinear |33, 44, 45], and high-energy be- 
havior [46, ^ which all become more evident in a helicity basis. The full color dependence 



is also helpful for exploring the general structure of infrared singularities [31 



In general, scattering amplitudes with massless gluon exchange possess infrared (soft 
and collinear) divergences. In dimensional regularization two-loop amplitudes generically 
contain poles in the dimensional regularization parameter e = (4 — Z))/2 up to These 
singularities have been predicted and organized into a compact form by Catani [31[ and 
cancel from final physical results. As is now standard, we use Catani's formula and color 
space notation to organize the helicity amplitudes into singular terms (which do contain 
finite terms in their series expansion in e), plus finite remainders. The precise form of the 
1/e poles was not predicted in ref. [^ for general processes at two loops. However, it 
is now apparent that these terms also have a universal structure depending only on the 
external legs, based on explicit calculation 14, 15, 16, 17, IS, 1£, ^ 

matching to resummations [|4 
collinear 



lii, 

and constraints on the functional form as momenta become 



A number of dimensional regularization variants are commonly employed for QCD loop 
calculations. The conventional dimensional regularization (CDR) scheme [^] is usually 
applied in calculations of amplitude interferences, such as in refs. [|5l|, |l^, |l^, 0, In 
the helicity approach, the two commonly used schemes are the 't Hooft-Veltman (HV) 
scheme 



and the four-dimensional helicity (FDH) scheme [53, 43 1. These schemes differ 
in the number of polarization components for unobserved gluons. The 't Hooft-Veltman 

2e virtual gluon states (as does the CDR scheme), whereas 

assigns 2 states. The FDH scheme 



(HV) scheme [52| contains 2 



the four-dimensional helicity (FDH) scheme [53, |4^ 
is related to dimensional reduction (DR) [^^, but is more compatible with the helicity 
method, because it allows two transverse dimensions in which to define helicity. A more 
detailed description of the differences between schemes, as well as a definition of the FDH 
scheme beyond one loop, may be found in ref. [43|. As in ref. [21| we find non-trivial scheme 
dependence in the finite remainders. 

The continuation of the Dirac algebra from four to D dimensions suffers from a variety 
of ambiguities. One well known ambiguity is in the continuation of 75 [52, 55[. Another 
lesser known ambiguity, which we also investigate in this paper, is tied to charge conjuga- 
tion. The appearance of an ambiguity is perhaps not surprising in hindsight given that like 
75 in D = 4 the charge conjugation matrix is a discrete product of 7 matrices, whose contin- 
uation to D dimensions is inherently delicate. In QCD it is a simple matter to sidestep this 
ambiguity by avoiding use of any charge conjugation identities that are valid only in four 
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dimensions. For M = 1 super- Yang-Mills theory, the situation is more complicated, be- 
cause of the Majorana nature of the gluino implying it is its own antiparticle. At tree level 
and one loop it turns out to affect the amplitudes only at 0(e) and is therefore irrelevant 
through NLO in the coupling constant. Moreover, it is does not appear even at two loops 
in the previously computed four-gluon or two-fermion two-gluon helicity amplitudes. This 
ambiguity first becomes relevant at NNLO when there are at least four fermions present. 
From Catani's infrared divergence formula, an ambiguity in 0{e) parts of a one-loop am- 
plitude necessarily feeds into an ambiguity in the two- loop amplitude starting at 0{l/e). 
Nevertheless, it is reassuring that this ambiguity turns out not to affect the finite remainder 
of the two-loop amplitudes as long as the different loop orders are computed consistently: 
We find that the entire ambiguity, including associated finite parts, may be absorbed into 



the Catani subtraction. Moreover, as noted in ref. |21|, when computing a physical process, 
at NNLO one does not need 0(e) contributions to the one-loop amplitudes because such 
terms always cancel. Thus the ambiguity completely cancels from physical quantities, as 
one may have anticipated from general considerations | |5^ . 

For theories without infrared divergences, one can straightforwardly add local coun- 
terterms to undo violations of supersymmetry by the regulator. For theories with infrared 
divergences the situation is more subtle. In any case, it is simpler and more elegant to use a 
regularization scheme which automatically preserves supersymmetry [^^. The supersym- 
metry preserving properties of the FDH scheme |53, 43 1 have been verified explicitly at two 



loops, for gluon-gluon and gluon-gluino scattering amplitudes |T^, 21 1. In this paper 
we study the supersymmetry identities satisfied by gluino-gluino scattering amplitudes. 
The ambiguities in the gluino-gluino amplitudes, however, causes a difficulty, since the 
precise value depends on the details of how the calculation was performed. However, after 
subtracting the singularities using the Catani formula, in the FDH scheme we explicitly 
show that the finite remainders, which are free of the ambiguities, do satisfy the expected 
supersymmetry identities. 

This paper is arranged as follows. In section § we present the helicity and color struc- 
ture of the four-quark QCD amplitudes. The color space structure of the divergent part 
of the amplitudes is reviewed in section ^. Section ^ contains the one-loop QCD ampli- 
tudes, which appear in Catani's formula for the divergent parts. The one-loop amplitudes 
are presented in a form valid through 0(e'^) since these are needed in Catani's formula 
for the two-loop divergences. The structure of the finite remainders after subtracting out 
the divergences using Catani's formula is given in section |5[ These finite remainders are 
tabulated in appendix The J\f = 1 supersymmetric amplitudes, the ambiguity in their 
value and the supersymmetry Ward identities, are discussed in section |6|. Some auxiliary 
functions needed for describing the scheme shifts are given in appendix ^. 

2. Helicity and color structure 

The three QCD processes considered in this paper are 

q(pi,^i) + q(P2, A2) Q(P3, A3) + Q(p4, A4) , (2.1) 
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q{pi,Xi) + Q(p2,A2) 



(2.2) 
(2.3) 



where we use a "standard" convention for the external momentum (pi) and hehcity labehng 
(Aj), i.e., particles 1 and 2 are taken to be incoming, while particles 3 and 4 are assumed 
outgoing. The identical quarks cases are easily obtained from these, as discussed below. 

We handle ultraviolet and infrared singularities using dimensional regularization. We 
consider a continuous set of schemes, labeled by a parameter characterizing the number 
of virtual gluon degrees of freedom circulating in loops. Specifically, when the trace of the 
Minkowski metric is encountered, we set 



(2.4) 



corresponding to 2(1 — eS^) gluon states in the loop. Setting 5^ = 1 corresponds to the 
HV scheme [^] , while setting = corresponds to the FDH scheme . 

The CDR and HV schemes have the same standard MS coupling constant, as(ju). The 
coupling in a general 5r scheme is related to this coupling at NNLO by [E^ 



a 



2tt 



+ 



12 



(1 - Sr) 



as(^) 
27r 



+ 



(2.5) 



Henceforth, for simplicity, we suppress the 6r index on as{fi). 

We work with ultraviolet renormalized amplitudes. The relation between the bare 
coupling a" and renormalized coupling as(/u), through two- loop order, is |31], 



/^O 



2e 



1 



as(/^) ^0 



27r 




(2.6) 



where /i is the renormalization scale, = exp[e(ln47r + '(/'(l))], and 7 = — V'(l) = 0.5772 . . . 
is Euler's constant. The first two coefficients appearing in the beta function for QCD, or 
more generally SU{N) gauge theory with Nf flavors of massless fundamental representation 
quarks, are scheme- independent, 

IICa - 4:TRNf _ nCl - (IOCa + QCF)TRNf 



bo 



6 6 

where Ca = N,Cf = {N^ - l)/(2iV), and Tr = 1/2. (Note that ref. 
/3o = W(2vr), /3i = 6i/(2^)2.) 

The perturbative expansion of the qq QQ amplitude is 

,(0) 



M 



qq- 



{as{fi),n;{p}) = 47ras(^) 



(2.7) 



uses the notation 



(2.^ 



+ 



27r ' 'ii^QQ 



(m;M) 



2tt 



M 



(2) 
qq^QQ 



(/i;M) + 0(a3(^)) 
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where {p}) is the L*^ loop contribution. The same type of expansion holds, 



of course, for the qQ — > qQ and qQ — > qQ amplitudes. Equation (^^) is equivalent to the 
following MS renormalization prescriptions at one and two loops. 



qq^QQ qq^QQ g qq- 



M 



(1) 

qq- 
(2) 

qq~*QQ 



S7^M 



(2)unren 
qq-*QQ 



qq- 



2e 



qq- 



We consider the following set of independent helicity configurations h 



(2.9) 
(2.10) 



h = 


1 


q{pi 


+) 


+ q{p2, -) - 


Q{P3, -) 


+ Q(P4,+ 


h = 


2 


q{Pi 


+) 


+ q{p2, -) - 


-^0(P3,+) 


+ Q{p4, - 


h = 


3 


q{pi 


+) 


+ 0(P2,+) 


g(P3,+) 


+ Q{P4,+ 


h = 


4 


q{pi 


+) 


+ 0(P2, -) 


9(P3,+) 


+ Q{P4, - 


h = 


5 


q{pi 


+) 


+ QiP2,+) 


9(P3,+) 


+ Q{P4,+ 


h = 


6 


q{pi 


+) 


+ QiP2, -) 


9(P3,+) 


+ QiP4, - 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 



Other configurations are simply related to these by symmetries. For example, the q{pi, —) 
amplitudes are obtained by parity (P), while the qQ — > qQ and qQ — > qQ amplitudes are 
related to qQ — > qQ and qQ qQ, respectively, by an overall charge conjugation (C). 
(Overall charge conjugation is unaffected by the ambiguity to be discussed in section |6^ .) 
In defining these helicity configurations we impose helicity conservation on the quark lines. 

The cases where both quark lines are identical, i.e. qq qq, qq qq and qq — > qq, 
can also be obtained from configurations ( 2.1lD -( ^ri6 ). For example, q{pi, +) + q{p2, — ) — > 
9(P3,~) + 9(P4,+) can be obtained by taking process ( 2.11| ) and subtracting the process 
( |2.14 ) with p3 and interchanged. The relative minus sign is due to the Fermi statistics. 
More generally we have for the qq qq process. 



M +- 



M +- 



-4 



where we use q'^ as a shorthand for q{pi, +) and so forth. 
The color decomposition of the amplitudes is given by 



/i = 1,...6, 



c=l 



where the elements in the color bases are 



(2.17) 



(2.18) 





= (5*4 (5*? 


TV[21 




h = 


1,2, 








= s\s\, 


h = 


3,4, 


TrW 


= (5*4 (5*? 


TV[2] 


= (5*? 5*4 


h = 


5,6 



(2.19) 
(2.20) 
(2.21) 
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In our previous calculations of the gg gg and qq — > gg processes [18, 21 1, traces of 
products of color matrices appeared. For consistency of notation, we have maintained here 
the "TV" notation for the color bases, even though there are no traces in the present case. 

The helicity-dependent, phase-containing factors Sh arise from evaluating the ampli- 
tudes in the spinor helicity formalism [E^]. They are. 



Si 
Sa 



'(4 2) 
'(3 2) 



S2 
55 



'(3 2) 

.(21) 
■(43) 



^3 
56 



.(21) 

'(43) 

.(41) 

'(2 3) 



(2.22) 



where '"'^^ 



The spinor inner products pTI , 33| are {ij) = ) and [i j] = {i' \J ), wnere are 

massless Weyl spinors of momentum ki, labeled with the sign of the helicity. They are 
anti-symmetric, with norm \{ij)\ = \ [ij] \ = ^/Si], where Sij = 2k_ 
prefactors Sh are all pure phases, i.e., 



kj. Notice that the 



\Sh\' 



1 



h = 1...6. 



r(0),[c] 



(2.23) 



The quantities M^"^'^''' depend only on the Mandelstam variables s = {pi + P2) > t = 
{pi —PiY and u = {pi — ps)^. At tree level in the color bases ( |2.19| ), ( |2.20| ) and ( |2.21| ) they 
are given by. 



M- 



(o),[i] 



M, 



(o),[i] 



(0),[2] 



M. 



(0),[1] 



N 



M, 



(o),[i] 



(o),[i] 



s 
s 



M. 



{o),[i] 



(o),[i] 



u 

t 

u 



/i= 1...6. 



(2.24) 



3. Infrared singularities 

We now briefly review the structure of the soft and collinear singularities of dimensionally 
regularized one- and two-loop QCD amplitudes, using Catani's color space notation pT|| . 
The finite remainders are given in section 4^, section 5^ and appendix |A[ 

The infrared divergences of renormalized one- and two-loop n-point amplitudes are 
given by [||, ||], 

|>1«(/x;M))h.. = I^'\e,f,;{p}) |M°^(^; {p}))h.s. + W)>aa , (3.1) 

\M^^Hf^;{p})h. = i^'\e,fi;{p}) |M'H/^;{p}))rs. 

+ I^il (e, /z; M) W (^; M)),s + [-^^'^'^'^(/u; M))k. , (3.2) 
where the "ket" notation \ {p}))rs. indicates that the L-loop amplitude is treated as 



a vector in color space. The components of this vector are given by the M^^^'^"^' in eq. ( 2.1^ ). 

The subscript as indicates that a quantity depends on the choice of regularization and 

renormalization scheme. The divergences of Ain^ are encoded in the color operator 

(2) (2) 
while those of A4n also involve the scheme-dependent operator 
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The operator which controls the one-loop singularity structure, is given by 



1 p-^i>W " 



i=i 



1 . jLi 

e2 + T? e 



(3.3) 



where Ajj = +1 if i and j are both incoming or outgoing partons, and Ajj = otherwise. 
The color charge Tj = \Tf^ is a vector with respect to the generator label a, and an SU{N) 
matrix with respect to the color indices of the outgoing parton i. For external fermions, 
the ratio 

7g_ ^ 3 
T? 2 ' 



(3.4) 



is independent of the representation. For quarks, = Cp = {N^ — 1)/{2N); for gluinos 
Tf = Ca = N. The two-loop operator I^H is ^| 

/i'^(e,/x;M) = -i/a)(e,;,;{p}) (^/a)(e,;.;{p}) + + 



+ 



+ l/L'l(e,/i;M), 
where the coefficient -fCaa depends on 5r and is given by [31, 15] 



(3.5) 



(2) 

The function -ff contains only single poles, and splits into two types of terms, 



67 ^2 



18 6 



6 9 



Ca - ^TnNf . 



(3.6) 



i/(2)( 



-eV(l) 



4er(l -e) 



2 \ 2e 



4Ff ) 1 + -H" 



(2) 



(3.7) 



where a standard analytic continuation is needed to bring s to the physical region. We find 
that the term proportional to the identity matrix in color space 1 is given by 



^(2) 



2 48 



23 o 245 



216 



CaCf + -6C3 + 



vr 



+ ( -^CaCf + \cl + \cFTRNf ] (1 - ,5^) 



vr^ _ 25 
12 ~ 54 



(2) 

This term survives the sum over colors, and the expression for Hq ' in the HV scheme 
{5r = 1) agrees, as expected, with previous color-summed results in the CDR scheme |13, 
|I|,||,|5|. 

The second term in H^"^^ (e) has exactly the same type of nontrivial color and kinematic 
dependence found in the gg ^ gg and qq gg helicity amplitudes |18, 21 1, namely 



CfTrNj 

(3.8) 



H 



(2) 



4 In — In In x 



Ti • To , To • T, 



(3.9) 



where again appropriate analytic continuation are required. For example, in the s-channel, 
ln((— s)/(— t)) — > Ins — ln(— t) — ivr. In refs. |18, 22| it was shown that for the gg —>■ gg 
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and qq gg amplitudes the structure of this term is independent of the hehcity configura- 
tion, and whether the external legs are quarks or gluons. Here once again, we find this to 

" (2) 

be the case. An ansatz generalizing H for an arbitrary number of external legs has been 

r(2) 



presented recently in ||4^. Because of the commutator structure, H vanishes when sand- 
wiched between tree amplitudes, after performing a sum over colors; hence it drops out of 
the color-summed interference of the two-loop amplitudes with the tree amplitudes 21 1. 
For each color basis we will have a different i"^^^ matrix. For the basis ( 2.1S| ) we have. 



where 



r(i 



s -u 



T-U 
2CpS-^(l 



U) 



(3.10) 



A* 



U 



1 3 
^ + 



2e 



(3.11) 



The corresponding operator for qQ —>■ qQ in the basis ( p.20D is obtained by changing S ^ U, 
T ^ S and U ^ T in eq. ( p.iop . Similarly, the operator for qQ —i- qQ in the basis ( 2.21 ) is 
obtained by exchanging S and U in ( p. 10 ). 

A typical partonic cross section requires an amplitude interference, summed over all 
external colors. Such interferences are evaluated in the color bases ( 2.19[ ), ( p. 20 ), ( 2.21| ) as 



j-{L,L') i^(i') V j^,j{LUc]* ^ ,MiL'),[c'] 

A1A2A3A4 \ A1A2A3A4I A1A2A3A4/ / J A1A2A3A4 A1A2A3A4 ' 



(3.12) 



c,c'=l 



where the symmetric matrix Q2c 




(3.13) 



4. One-loop QCD amplitudes 

The one-loop amplitudes for qq — > QQ were first evaluated through 0{e^) as an interference 



with the tree amplitude in the CDR scheme |51]. Later they were evaluated as helicity 



amplitudes in the HV, FDH and DR schemes [30|. 

Because I^^^ contains terms of order 1/e^, the I^^^\Ai^^^)^s. term in the infrared de- 
composition ( ^ ) of the two-loop qq — QQ amplitudes requires the series expansion of 
the one-loop amplitudes through O(e^). In section 4.1 we present the all-order results 
in the color bases ( 2.19| ), ( |2.20 ), (|2.21| ), with the normalizations implicit in eq. ( |2.8[ ), in 



terms of integral functions whose series expansions have been evaluated to the required 



order [11, 12]. 



In ref. |18| it was shown that the 0{e) terms in one-loop amplitudes such as A4 . 
are not required for the construction of a numerical NNLO program, once the divergences 



have been subtracted from M 



(2) 



using Catani's formula. Thus we need only present 



explicit formula for the 0{€^) finite remainders Ai 



,(l)fin 
qq^QQ 



of the one-loop amplitudes, after 



8 



ultraviolet renormalization ( p.9| ) and subtraction of infrared divergences ( |3.1| ). We do this 
in section 4.2, for a general Sn scheme. The formula (5.10) for converting the two-loop 
finite remainders A^^^J^^„ from one scheme to another are most compactly presented in 
terms of the J/j-dependent parts of the one-loop amplitudes at order e; the explicit values 
of these quantities are collected in appendix 



4.1 One-loop all orders in e QCD amplitudes 

We now present the one- loop qq — > QQ amplitudes in the color bases (p.l9|) , (|2.20| ), (|2.21j ), 
with the normalizations implicit in eq. (p.§|), in a form valid to all orders in e. 

At one loop the crossing properties of the amplitudes are relatively simple, so we 
present the explicit values of the helicity amplitudes for the process qq QQ. The process 
qQ qQ may be obtained from these by crossing the initial antiquark and final quark 
into the final and initial states respectively, and the process qQ — > qQ may be obtained by 
crossing antiquark 2 into the final state, and antiquark 3 into the initial state, 

r(l),W, 



M, 
M 



4 



{s,t,u) = Ml'>'^''\u,s,t), c 
{s,t,u) = M[^^'^''\u,t,s) , c 
{s,t,u) = M^^^'^''\u,t,s) , c 



1,2, 
1,2, 
1,2, 
1,2, 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where M^^^''"' is defined in eq. ( |2T8| ) with the color bases (|1|), ( ^3oD and ( ^) using the 
helicity configurations h defined in eqs. (2.11)-(2.16). After crossing, appropriate analytic 
continuations are required to bring each function into the physical region. 
A compact representation for the unrenormalized amplitudes is, 



M- 



{1),[2] 



M. 



{1),[2] 



{s,t,u) 
(s, t, u) 
{s,t,u) 



A+{s, u, t) + —Ai{s, t,u)+ (N - — ] A2{s, t, u) 



N 
{s, u, t) 



^ + 7^ )Ai(s,t,u) + -^A2{s,t,u 



iV2- 



A_{s, t, u) + —Asis, t,u) + [N - — ] A4{s, t, u 



1 

N 



N 



A^{s,t,u) - 1 + 



N 



1 

iV2 



1 



A3{s,t,u) + j^A4{s,t,u) 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



where, 

A±{s,t,u) 

Ai{s,t,u) 

A2is,t,u) 



2e- 1 



e(l 



2(e-l) V2e-3 
1 

~ N 



2 2e 

N + 



1 



2e-3 



N 
Tri(4)(^ 



+ 



2(2e - 1) 



u 



2(2e 



—tSr{1 + e6R)u H 

1 s 



Tri(^)(u) 



y(5i?(l + e6R)s 
sTri(^)(s) 



e'6R -{2e-l)- ]u 



Box('^)(s,u), 



2(2e - 1) 



^TViW(i) + 



(4.9) 



(4.10) 
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^4(5, t, u) 



+ 



il + e'6ji)s- [e'5n + {l-2e)- ]u 



Box(^)(s,i) 



(4.11) 



-(e - 1)6r{1 + eSR)s + e^6Rt - e- + {1 - e) 
2 s 



+ 



'-6R{l + e6R) + {l-2e)i-^ 



u 



2e- 1 



2s 



2(2e - 1) 

^3 



2(2e - 1) 

.3 



Sit 



t 



2i 



Tri(^)(t) + 



+ 



^dls - e^6Rt + {{I - e)u + et)- 



Box(6)(s,t). 



(4.12) 



(4.13) 



Here Tri*-^^ (s) is the scalar triangle integral in 4 — 2e dimensions with one external massive 
leg, and Box*-^^(s,t) is the all-massless scalar box integral in 6 — 2e dimensions. The 
expansion of these integrals to Ole^) in the various kinematic channels is given, for example, 
in refs. |jl2|, |l6[. The renormalized amplitudes are obtained by subtracting boM^^^'^'^^ /e from 
each of eqs. (E^-(H?8|), where c, h correspond to the amplitude under consideration. 



4.2 Finite remainders 

We now give the finite remainders of the one-loop qq — > QQ, qQ qQ and qQ qQ 



amplitudes at 0{e^), defined by Ai^g^^gg and -Mqg^^gq in eq. ( ^ ) and color decomposed 



(l)fin 



into M; 



(l),[c]fin 



m eq. (|2l8D. We write. 



(l),[l]fin 



60 I In 



ivr + 



N 



^) 



M, 



(l),[2]fin 



60 I In 



iir + 



N 
3" 



1 
2N 



(I-Sr) 



^(o),[i] ^ 



, [2] [2] 1 ,[2] 

+ ^^h + ^2 % 



(4.14) 



(4.15) 



For the h = 1 helicity amplitude, the independent remainder functions a, b, d, h, j 
and k are 
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X 



1 



—y+[--y]X + -[y + — ]X' + iTT 



X 1 f x"^ . ^ 



2\y 



b['^ = (2y - x)X 



l/x' 



+ y)x'^ + y(Y^ -3Y + + 



+ iir 
5 



2\y 



+ y]X + 2yY 



(4.16) 

(4.17) 
(4.18) 
(4.19) 
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9^' 



1 fx' 



+ y]X'' + y[-Y 



where 



— ZTT 



X 



4 V 2/ 

^K7 + ^)^ + ^^^ 



4 - 



s s 
For /t = 2 the functions are 

13 



X = ln{-x) , Y = In(-y) . 



a. 



[1] _ 



X' -3X + 



+ i-KX X 



^ fy 



x{A + ^X - X^) + - [ — + x]Y^ + {y - 2x)Y - 



— ZTT 

5 

-X . 



2\x 

l + 2xX- [x + —\Y 

X 



— — X ]Y — iir 



^-x + -{x + ^]Y 



X 



9" 

4 V X 

For h = 3 the functions are 



x)y + a;( -X' 



4--X)+ 



1. 



1. 



13\ 1 



18; y 



ivr y 



3/!/ 



TT 



(1 + x^) - - Xy + -X^ + 4 - (2 - x)X + - {x' - 1)Y^ 



- + Y- 

y 



-i7r(2y + 3)- 



y 



4' 



[21 



,•[2] 
^3 



= -y- 

3 

(l + x^ 

1 



5\ 1 



vr 



y 3y 

ixy-i(7r2 + y2 + x2' 

2 4^ 



X 



+ 1-- x + -y- 



13 
18 



-y + ivr— 
2 6y 



1, 



-y 



TT 



(1 + x')i2XY - TT^ - X") - 4 + ( 1 - ^ )X + ^(1 - x2)y2 



1 1. 



(4.20) 

(4.21) 
(4.22) 

(4.23) 

(4.24) 
(4.25) 
(4.26) 

(4.27) 

(4.28) 
(4.29) 

(4.30) 
(4.31) 

(4.32) 
(4.33) 



+ ivr y + 



3\ 1 



2Jy 



(4.34) 
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For h = A the functions are 



6W 



r /i3 
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ITT 
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X llx 
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ITT- 
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+ iir 
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For h = 5 the functions are 
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18 ~ 
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1.. ■ 11 

2 6y 
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y 
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3 9J y 3y 



1. 



/^r^ = -^y 



.[2] 

J5 



1. 



J'-Y-'-y--^^ 
\3 9jy 3y 

1 



-(1 + x')iX' + 2XY) - 4+ (- -l)X--(l- x')Y 
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For h = 6 the functions are 
-2 3 „ 1 



a, 



[1] 

6 



2 2 



-Y + -(X^ + y'^)-XY + — 
3 2^ ^ 18 



x(4 - TT^ + 3X - + 2Xy) + 



1 n 



2 Vx 



X y 



X llx 

ITT , 

y 6y 

i + r + 
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(4.35) 

(4.36) 
(4.37) 
(4.38) 
(4.39) 

(4.40) 



(4.41) 

(4.42) 
(4.43) 
(4.44) 
(4.45) 

(4.46) 

(4.47) 
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ITT- 
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^( 
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(4.48) 
(4.49) 

(4.50) 
(4.51) 

(4.52) 



For the HV scheme {b^ = 1), the results (4.14)-(4.52) for the finite remainders of the 
one- loop helicity amplitudes are in agreement with those of ref. [p^ j. 

5. Two-loop QCD amplitudes and finite remainders 
5.1 Construction of amplitudes 

We generated the Feynman graphs for qq — > QQ using QGRAF pl|, from which a MAPLE 



program was constructed to evaluate each graph. We employed the integral reduction 



algorithms developed for the all-massless four-point topologies p9| , [30| , 28, 27, in order 
to reduce the loop integrals to a basis of master integrals. To put the integrands into a 
form suitable for applying the general reduction algorithms, spinor strings were converted 
to traces over 7 matrices, by multiplying and dividing by appropriate spinor inner products 
constructed from the external momenta. To illustrate the method consider the diagram 
depicted in fig. |l[ The numerator of the integrand is 



We multiply and divide this by 

(l-|4+)(3+|2-) = (14)[32] - ^^^^^ 
so that the numerator can be rewritten as, 



(23) ' 



(5.1) 



(5.2) 



^^^^^-(2-|7^'fe7'^^27^^i7/.|l")(l"|4+)(4+|7p^47.|3+)(3^ 



(14) t 



^(2|P+7^47'^l^27''^i7^^-|l)(l|4)(4|P_7pl^47.^'+|3)(3|2) 



spins 



Mi 

(14) t 

MiTr(^2^+7^/^37"^27'^i7M^-^i^4^->^47.^+|^3 
(14) t V^2 ' 



If 4 



P_7p|^47,,P+|i, 



(5.3) 



13 



(3+|2-> 




(l-|4+> 



(l-\fA2-) IJ 




(b) 



h (3+|^J4+) 



Figure 1: An example of a two- loop diagram showing two different ways of "closing" the fermion 
lines. The dashed lines represent the spinor inner products inserted in order to form a trace. 



where pi is the momentum of external leg i and we have introduced the helicity projectors 
P-f and P_, given by 

^'+ = ^(1 + 75), P- = ^(l-75), (5.4) 
in order to reduce the sum over spins to the original helicity configuration. 

To evaluate the projectors we use the 't Hooft-Veltman prescription [^] for 75, 

{75,7!'^} = 0, [75,7(-'^)]=0, (5.5) 



where the notation "(4)" and "(— 2e)" is used to indicate whether the Lorentz index fi lies 
the four-dimensional or (— 2e)-dimensional subspaces, respectively. (Another prescription is 
to take 75 to anti-commute |5^ with all components, but this has the unwanted side-effect 
of ruining the smooth connection of the 't Hooft-Veltman and FDH scheme as a function 
(5r.) Using eq. ( ^ ) we move one of the projectors in the last line of eq. ( ^ ) until it hits 
the other, giving, 



Ml 

(14) t 



Tr 



,7' 



■Tr(^2^^^37'^^27"^i7M^i^47i-'^^ri'^^'^^^+^: 



+ 
+ 



(5.6) 



In the last equation, the 75 part of P-f will produce terms containing a Levi-Civita tensor. 
Upon integration this tensor can appear contracted only with external momenta, and since 
we have only three independent ones, these terms must vanish. Therefore the 75 term from 



the P+ in eq. (5^) can be dropped here. 

After expanding the traces over Dirac matrices using standard formulae, we obtain dot 
products of momenta (external and/or internal) in our expressions, making possible the 
application of the reduction algorithms previously mentioned. In evaluating the trace we 



14 



take the four-dimensional and (— 2e)-dimensional subspaces to be orthogonal following the 



discussion of ref. |43]. The rules for evaluating integrands containing dot products of (— 2e) 
dimensional components of momenta may be found in ref. 

Alternatively, we could have also chosen to "close" the fermion lines in a different way, 
multiplying and dividing, for instance, by 

(l-|^3|2-)(3+||)j4+) = (14)[32]n. (5.7) 

This is depicted in figure |l|(b). However, this would have given a product of two traces 
instead of just one, which would produce some terms containing a product of two Levi- 
Civita tensors, which would then need to be evaluated since they would not vanish under 
integration. 

After all the tensor loop integrals in the amplitudes have been reduced to a linear 
combination of master integrals, the next step is to expand the master integrals in a 



Laurent series in e, beginning at order 1/e , using results from refs. |25, 26, 2S 



The results can be expressed solely in terms of ordinary polylogarithms [32, 63| 

'■^ dt 



Unix) 



Li2(x) 



oo i 

Ex 
An 



=1 







t 



Li„_i(t) 



(5.8) 



(5.9) 



with n 



2, 3, 4. The analytic properties of the non-planar double box integrals appearing 

|; there is no Euclidean region in any of 



in the amplitudes are somewhat intricate |11, 
the three kinematic channels, s, t or u. We do not attempt to give a crossing-symmetric 
representation, but instead quote all our results in the physical s-channel (s > 0; t, u < 0) 
for the qq — > QQ, qQ — > qQ and qQ qQ kinematics, eqs. (2J.), ( |2.2| ) and (p.3[). 



5.2 Finite remainders 

The two-loop finite remainders are defined in eq. (^^) and are color decomposed into 



{2),[c]fin . 



h 



in eq. ( 2.18| ). Their dependence on the renormalization scale fx, the number of 



colors A^, the number of fermion flavors A^, and scheme label may be extracted as 



M. 



{2),[c]fin 
h 



where, 



-hi (ln(s/^2) _ .^)2 _ j,^ (in(s/;,2) _ - (Ca + QCf)Ca{1 - Snf + 



36 



+ {4i?g + ^{Ca + 9Cf)Qo + 6oQS''^(ln(s//x') 



in) 



1 



+ 



JCA-6CF)QoiTr){l-5R) 
-2boiHs/^l^)-^n) + Q^^'^\l - Sr) 



Mmc] + 

^(l),Wfin ^ 



+ QoM«'W^'^-(1-5h)+41 



p, 



[1] 



Nf 



p, 



[2] 



(5.10) 



(5.11) 



15 



The //-dependence is a consequence of renormalization group invariance. 

The tree and one-loop functions, Mj^^'^'^^ and Mj^^'^'^^^^, are given in eq. ( |2.24| ) and 
eqs. (|4ll and ( p^ ), while 60 and 61 are given in eq. ( |2.7D . The following combinations 



of color constants also appear in eq. (|5.1[l|) , 

Qo = ^Ca -Cf + ^TnNf , (5.12) 
q[^i) = _^Ca + Cf, (5.13) 



The constants Qo and Rq also appeared in the two-loop finite remainders for the qq — > gg 
process, and Q^'^^ is just twice the constant Q^^^ appearing in that case. The quantities 
jy^WJc] fi'Jfl ^j^^ (5^-dependent parts of the O(e^) coefficients of the one-loop amplitude 
remainders, after subtracting the poles using eq. ( |3lD . The explicit values of M^^^'^"^ "'"^^ are 
tabulated in appendix ^. The coefficient functions A, B, C, D, E, F, G, H, I, J, K depend 
only on the Mandelstam variables. In appendix we give the explicit forms for the 
independent finite remainder functions appearing in eq. ( ^.11 ). 



We have compared our results for the independent two-loop finite remainder functions 
j^^2),[c]fin ^-^.j^ corresponding results obtained recently by Glover ref. [23|. Our results 
agree with the corrected version of ref. [23|, once a slightly different definition of H^'^\e) in 
eq. (|3.7D is accounted for. (Instead of dressing the 1/e pole with (/x^/(— s))^*^, as we do in 
eq. in ref. |2| it is dressed with {n'^/{-s)f^ + {n^/{-t)f' - {n^/{-u)f'). We have 

also checked that the interference of the tree and two-loop helicity amplitudes, summed 
over helicities and colors, reproduces the results given in ref. [p^ j. 

6. Amplitudes in pure = 1 super- Yang-Mills theory 

Supersymmetric gauge theories have a wide range of applications both for phenomenological 
and theoretical purposes. Here we present the amplitudes in = 1 pure super- Yang- 
Mills theory, obtained from the QCD ones by modifying the fermions to be in the adjoint 
representation and by altering their multiplicity to correspond to a single Majorana fermion 
circulating in the loops. (A Majorana fermion counts as half of a Dirac fermion.) The 
fermion in the M = 1 pure super- Yang-Mills theory is the gluino superpartner of the 
gluon. Besides the inherent interest in supersymmetric theories, a useful consequence is 
that supersymmetry imposes a set of powerful constraints on the amplitudes which can be 
used to verify their correctness. The supersymmetry identities have been applied previously 
to the same one-loop amplitudes discussed here [30|, but only through 0{e^), as needed 



in an NLO calculation. In this section we discuss the supersymmetry identities for gluino- 
gluino scattering up to two loops. 

6.1 Supersymmetry Ward Identities 



In refs. ||l^, ^ it was shown that using the FDH scheme {6r = 0) the following identities 
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(a) (b) (c) 

Figure 2: Three contributions for Majorana fermions. The assignment of fermion arrows is incon- 
sistent since there is no distinction between particles and antiparticles. 

are satisfied through 0{e^) at two loops, 

M''''''''igt,92,9t,9t) = 0, (6.1) 
A^'^'^ (5^52^53^,54+) = 0, (6.2) 

-M'^'^(5^52-,53,5|) = ||-M'^'^(5^52-,53-,54+), (6.3) 

where g and g denote a gluon and gluino, and the superscripts denote helicities. There is 
also an identity relating the four-gluino amplitudes to the four-gluon ones, 

A^^^^^(5^52-,53-,54+) = -^^M'^^^igt,g,,gs,9t) 

= -^^M^^'^rgt,9^,9s,9t) (6-4) 
which we discuss here. In ref. |60| this was shown to be valid at one loop through 0(6*^), 



as required in an NLO calculation. 

6.2 Ambiguities in D-dimensional Dirac algebra 

Before presenting the results, we first describe ambiguities that affect the amplitudes. 
The methods we use rely on continuing the Dirac algebra appearing in the spinor inner 
products away from four dimensions. This continuation is of course not unique, with a 
variety of schemes such as the 't Hooft-Veltman or FDH scheme for doing so. However, 
even within each of these schemes there can be further ambiguities. In particular, the 
charge conjugation properties of inner products of helicity states are not well defined as 
one moves away from four dimensions. In four dimensions we have the identity, 

{Plhf^l ■ ■ ■ 7M2n + l Ipf) = (P2 l7M2n + l " " " 7^1 \P2) , (6-5) 

which is useful when performing Fierz rearrangements of the spinor products. It can also 
be used to construct simpler Dirac traces to evaluate. But can we use this identity when 
we continue the Lorentz indices of the Dirac matrices away from four dimensions? 

For the case of Dirac fermions such as the quarks of QCD, we can avoid answering 
this question by imposing the rule that eq. (|6.5D should not be used in any manipulations. 
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Indeed, in performing the QCD calculations described in the previous section, we did not 
use this equation. However, with Majorana fermions, as appear in the N = 1 theory, there 
is no distinction between particles and antiparticles. In particular, one cannot consistently 
enforce a global requirement that a given fermion be treated as a particle and no global 
assignment of fermion arrows is consistent for all diagrams |6^. For example, in figure ^ 
if we chose legs 1 and 4 to be particles and legs 2 and 3 to be antiparticles, this would be 
inconsistent with diagram figure |2|(c) which is also a perfectly valid diagram for Majorana 
fermions. 

Consider, for example, the box diagrams of figure ^. With Majorana fermions it is not 
clear which of the two diagrams should be used. The two diagrams differ by an application 
of the identity ( |6.5D and have identical values in four dimensions, up to the overall sign. 
In order to calculate these diagrams we can multiply and divide the integrand of diagram 
(a) by (l~|4+)(3+|2") and the integrand of diagram (b) by {l~\f^\'i~){A~\f^\2~). After 
performing the Dirac algebra, it is straightforward to show that the integrands of these 
two diagrams differ by a term proportional to A^, where A is the (— 2e)-dimensional part 
of the loop momentum. Since the one-loop box diagram with a single A^ in the numerator 
is of 0{e), this ambiguity is irrelevant at NLO. However, at two loops this ambiguity does 
enter, as one might expect by observing that 0{e) terms at one loop contribute to 0{l/e) 
terms at two loops in Catani's formula ( |3.2D . With only a single or no fermion pair the 
results are unaffected by applying eq. ( |6.5p . This ambiguity therefore does not appear in 
the gggg or gggg amplitudes. 

A related issue is the lack of helicity conservation at a fermion vertex with the 't Hooft- 



Veltman prescription |52| for 75 (5.5). Since the (— 2e)-dimensional part of 7-matrices 
commute, helicity violating terms of the form (a~|7^~^'^^|6^) = (a|P+7(~^'^)p+|6) do in 
general contribute. At tree level the helicity violating contributions vanish in the FDH 
scheme. However, they do contribute at loop level because 0{e) terms are generated which 
interfere with the divergences. This is a reflection of the well known violation of the chiral 
Ward identity when using the 't Hooft-Veltman 75 prescription. Because of their connection 
to the divergences, it is not surprising that the fermion helicity violating terms can all be 
absorbed into Catani's formula for two-loop divergences, as confirmed by our calculations, 
and hence do not affect physical quantities. It is therefore perfectly consistent to drop 
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these contributions (as we have done in QCD), if this is done systematically throughout 
the calculation. In = 1 super- Yang-Mills theory, although it is again consistent to drop 
these contributions, they do affect super symmetry Ward identities. 

One can choose various prescriptions for resolving these ambiguities, but it is not a pri- 
ori clear which will allow the amplitudes to satisfy the supersymmetry identities. Whether 
or not one keeps fermion helicity violating contributions and the choices for assigning 
fermion arrows alters the value of the full amplitudes. Effectively, these ambiguities mean 
that the supersymmetry Ward identities will not hold unless an additional prescription is 
imposed. As we discuss below it is possible to choose the ambiguous terms so that the 
supersymmetry identities are satisfied. 

However, a more straightforward way to deal with these ambiguities is to use the fact 
that they are all tied to the divergences controlled by Catani's formula. Our explicit cal- 
culations verify that the finite remainders after subtracting the divergences using Catani's 
formula are completely free of these ambiguities, as long as the tree, one-loop and two-loop 
amplitudes are treated uniformly with the same prescriptions. That is, we obtain the same 
finite remainders whether or not we apply eq. ( |6.5D to one of the fermion lines, or whether 
we keep or drop contributions that violate fermion helicity conservation. Moreover, as we 
discuss below, in the FDH scheme the finite remainders satisfy the expected supersymmetry 
identities without any additional prescriptions. 

6.3 Color and infrared structure 

Because the four gluinos are identical particles, and are Majorana, if we enforce helicity 
conservation, there are only two independent processes, 

h = l: g{pi,+) + g{p2,-) ^ g{p3,-) + g{p4,+) , (6.6) 
h = 2: g{pi,+)+g{p2,+) ^ g{p3,+)+giPi,+)- (6.7) 

We ignore the helicity violating processes, e.g., 

h = Z: g{pi,+) + g{p2,-) ^ g{p3,+) + g{p4,+) , (6.8) 
h = 4: g{pi,-)+g{p2,-) ^ g{P3,+)+9{P4,+), (6.9) 

since, as mentioned above, and confirmed by explicit computation, they can be absorbed 
entirely into Catani's formula for divergences. 

Since the gluinos are in the adjoint representation we use the same color basis as used 
for the four-gluon helicity amplitudes fl^ 



Mi"-^ = ShxY, Trl'^l X Mf ^'['^1 , (6.10) 

c=l 

where 
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A reflection identity implies that the c = 4, 5, 6 coefficients are equal to the c = 3,2,1 
coefficients (respectively), so there are really only six different coefficients for each h, namely 



~AL),[c] 



h 



, c = 1, 2, 3, 7, 8, 9. The corresponding spinor factors are 



— 



.(31) 



So — 5*5 



:(21) 



(6.12) 



(42)' "(43) 

The infrared divergence structure is similar to that of gluon-gluon and gluino-gluon scat- 
tering amplitudes 21 1. For the case of = 1 pure super- Yang-Mills theory, in the 
basis (KW) the matrix is P, pi, l2l 



^ (^) = -f(T^b + iv^, 

/7V(S + T) 

7V(S + U) 

7V(T + U) 

iv(T + u) 







(6.13) 



(S - U) (S - T) 


V (T-U) 
















iV(S + U) 


(S-T) (S-U) 



(T-U) (S - U) 

(U-T) (S-T) 

(T - S) (U - S) 

(T - S) (U - S) 

(U-T) (S-T) 

iV(S + T) (T-U) (S-U) 



2NS 

(U-T) (U-S) (U-S) (U-T) 27VU 

(T-S) (T-S) (T-U) 2m J 



where S, T and U are defined in eq. ( p.ll| ). For = 1 super- Yang- Mills theory the ffi'st 
two coefficients of the /3-function are 



^0 = -^Ca , 



(6.14) 



The /'•^^ operator for 

f(2) 



super- Yang-Mills theory is 
^ e+^^Wr(l-2e) 



+ 



r(2) / r ^^ 



where 



/ 4 \ 



■(2)^ 



(6.15) 

(6.16) 
(6.17) 
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and 



/ 7_7'(2)\SYM 



(2)nSYM 



FDH 



— + — 

2 16 



r2 



(6.18) 



The equality of (^rf^)|^^ and {H^'> 

)fdh ^ consequence of supersymmetry. Equa- 



r(2)^SYM 



tions ( |6.16| )-( |6.18 ) are obtained from the QCD formulas, eqs. (3.7)-(3.8), by the replace- 
ments 5ji — > 0, Cp — > Ca and T^Np — > CaI'^ for converting to a single adjoint fermion in 
the FDH scheme. The operator H defined in eq. ( p.9[ ) does not explicitly depend on the 
fermion representation. 

The tree amplitudes in this color basis in the FDH scheme are, 
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and are free of the ambiguities described above. 

6.4 One-loop amplitudes in pure M = 1 super- Yang-Mills theory 

We now present the results for one-loop four-gluino scattering in a format valid to all orders 
in e. However, because of the ambiguity discussed above the answer depends on the precise 
steps used to calculate the expression. 

We can, for example, form the h = 1 gluino amplitude from the h = 3 (with traces 
formed by multiplying and dividing by [12] (43)) amplitude ( |2.13| ) with adjoint represen- 
tation quarks and Nj = 1/2 summed with its 3 <-> 4 interchange. This gives, for the 
unrenormalized coefficients of the three color structures Tr'^l , Tr''^! and Tr^^l , 
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(6.20) 
, (6.21) 
, (6.22) 



where we have taken the FDH scheme {Sr = 0). The renormalization is trivially performed 
by subtracting the quantity h^M^^'^''^ / e from each, where 6o is given in eq. ( |6.14 ). The 
ambiguities are all proportional to eBox*^^^ and hence are of 0{e). 

Other choices are also possible. For example, if include a fermion helicity violating term 
q^q^ — > Q'^Q^ (forming traces by multiplying and dividing by (4~||^3|l~)(3~|/^i|2~)) in 
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constructing the h = 2 gluino amplitude, this shifts eq. ( |6.21 ) by 



5m\^'^'^^^ = ^iVneBox(^)(n,s) , (6.23) 

leaving the coefficients of the first and third color structures unshifted. Similarly, adding 
in the fermion helicity violating contribution q^q^ — > Q^Q~^ shifts eq. ( 6.22| ) by 



For any of the above contributions, we can swap one of the particles with its antiparticle be- 
fore converting them to gluinos, which again alters their values by 0(e) terms proportional 
to eBox(^). 

The remaining color coefficients, up to 0{e) ambiguous terms, are given in terms of 
the previous ones by, 

M«'[^-^] t, u) = M«'[^+^] (s, t, n) i = 0, 1, 2, , 



(s, t, u) = M^'l^l (s, t, u) = M^^^'f^l (s, t, u) , (6.25) 

The h = 2 amplitudes are related to the h = 1 through the following relations, 

M^^^'W(s,t,n) = -Mf)'[2l(t,u,s), 
Mi^)'['l(s,t,u) = -M['^'^'\t,u,s), 
M«'['l(s,t,u) = -Mf^'''l(t,n,s), 

M2^^)'[^] (s, t, u) = -Mf (t, n, s) . (6.26) 
To check the supersymmetry identity (|6.4|) we compare these results to the correspond- 



ing ones for gluon-gluon and gluino-gluon scattering given in refs. We see that 

eqs. (|630|) , ( |63l|) and (|]2|) match eqs. (5.28), (5.29) and (5.30) in ref. |2| respectively, 
up to an overall normalization factor of —st/u^ (which is similar to the same factor in eq. 
(5.39) of ref . |pl[| ), except for terms proportional to eBox^^^ in M2^^''^^ and M2^^'^^\ (Note 
that in ref. [^] the amplitudes include the renormalization subtraction, not explicitly in- 
cluded here.) Similarly it matches the corresponding equations in section 3.1 of ref. [jl8| for 
gluon-gluon scattering, again up to the terms proportional to eBox^^^. Since Box^^^(s,ii) 
is finite, the one- loop terms that violate the supersymmetry Ward identity are of 0(e). 
These terms are precisely the ambiguous terms described above. We can alter such terms 
by modifying the prescription. If, for example, we include in the shifts ( |6.23 ) and ( 6.2'j| ) 



arising from fermion helicity violating terms we find that the supersymmetry identities are 
then satisfied to all orders in e. Of course, this is not completely satisfactory because of 
the ad hoc nature of such choices. 

6.5 Two-loop amplitudes in = 1 super- Yang-Mills theory 

The dependence of the identity ( |6.4| ) beyond 0{e^) on ambiguous terms at one loop foretells 
a similar dependence at two loops except this time at 0(l/e). However, after subtracting 
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these divergences, via Catani's formula, we find that eq. ( p.4| ) is indeed satisfied for the finite 
terms which are independent of the ambiguous terms and any prescriptions chosen for fixing 
them. However, some attention is required to ensure that precisely the same prescriptions 
are applied to tree level and one loop as applied to two loops. Otherwise, there would be 
a mismatch between the Catani subtraction terms and the two-loop divergences leaving 
behind dependence on the ambiguous terms. 
Our results in the FDH scheme are 
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c = 7,8,9. (6.28) 



Because the adjoint color indices of the gluino fields are identical to those of gluons, and only 
structure constants /"^^^ appear in the two-loop Feynman diagrams, the color coefficients 
for two-gluino two-gluon scattering obey the same group theory relations identified for 
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We have verified that the finite remainder functions match those of the pure gluon 



case, t.c 
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(6.33) 
(6.34) 

where X G {A, B, G}. The functions x^^j}^ and x^^_^^ for gg gg scattering in pure 
M = 1 super- Yang- Mills theory are given in ref. [^]. (Note that in that reference an all 
outgoing definition of helicity is used.) The relations ( |6.33| )-( |6.34| ) are precisely equivalent 
to the content of the supersymmetry Ward identities E^] , after removing overall factors and 



the divergent terms, and extracting the N and fi dependence. As discussed in ref. |21| this 
also matches the corresponding functions for the gg gg amplitudes. This demonstrates 



that the supersymmetry identity (6^) holds for the finite remainders at two loops. In the 
HV scheme, as expected, it does not hold because of the mismatch of fermionic and bosonic 
states. 
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7. Conclusions 



In this paper we presented the two-loop helicity amphtudes for quark-quark and anti-quark- 
quark scattering in QCD and gluino-gluino scattering in M = 1 super- Yang-Mihs theory. 
These amphtudes retain the fuh information on external colors and helicities. We verified 
that the interference of our two-loop amplitudes with the tree- level amplitudes, summed 
over all external colors and helicities, and converted to the CDR scheme, are in complete 



agreement with the results of ref. |13]. We also found complete agreement with recently 
published results by Glover [23|, after correction of minor errors in the original version of 
that article. 

We also discuss ambiguities in defining the amplitudes, related to the continutation of 
75 or charge conjugation identities away from four dimensions. As confirmed by our explicit 
calculations these ambiguities drop out of final physical results, since the ambiguous contri- 
butions can be absorbed into infrared singularities that cancel from physical quantities. In 
particular, there is no physical content to contributions that violate helicity conservation 
of a given massless fermion. However, some attention needs to be paid so that the differ- 
ent loop orders are computed consistently with the same set of prescriptions throughout. 
Otherwise, the cancellation of the ambiguous unphysical terms would not be complete. 

In ref. it was shown that the supersymmetry Ward identities relating the gg — > 
gg to the gg gg processes are satisfied in the FDH scheme at two loops, for both 
infrared divergent and finite parts. For the gluino-gluino scattering amplitudes discussed 
in this paper, unless the ambiguities entering in the amplitude are carefully adjusted, the 
supersymmetry identities will not hold starting at 0(1 /e). In any case, after subtracting the 



divergences using Catani's formula |31], the finite remainders are free of these ambiguities 
and satisfy the expected supersymmetry identities in the FDH scheme. 

So far, the new 2^2 amplitudes amplitudes have been implemented in a handful of 
new phenomenological studies []T7| , 35 1. Once general algorithms for dealing with infrared 



divergent phase space integrations at next-to- next-leading-order are completed fG^], many 
more phenomenological applications will follow. These applications would include the 



implementation of the two-loop four-quark amplitudes of this paper, or those of refs. |13 
|23[| , as ingredients in a numerical program for computing dijet production cross sections at 
hadron colliders at NNLO in QCD. When this task is accomplished, the intrinsic precision 
on the QCD predictions should reach the few percent level. 
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A. Finite remainder functions for QCD 

In this appendix, we present the explicit forms for the independent finite remainder func- 
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tions for the processes qq — > QQ, qQ qQ and qQ — > qQ in QCD, which appear in 
eq. ( pll) . 

For the heUcity h = 1 configuration in eq. (2.11) and color factor Tr^^^ in eq. ( 2.1S| ), 
the finite remainder functions are: 
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+ (^(2 + x)X + ^ + 2Y^ Li3(-x) - ( 



1 1 



24 12 



^ _ 1 
12 ~ 6 



xY]X^ 



+ y^x - 2y + - X + -TT^ + -tt^x Li2(-x) + 
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+ 
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For /i = 2 in eq. (2.12) and color factor Tr^^' in eq. (2.19): 
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[1] 
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18 



3Li4( -- 

y 



(3 - x)Li4(-x) + -xX2Li2(-x) + 
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11 
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30 
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+ 
+ 
+ 
+ 

+ 

+ 
+ 
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^X + 3XU2i-x) 
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— (5a; - 7)X^ + 12Li4 (^-- j + (4a; + 12)Li4(-x) + 2xX^Li2(-x) + 
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For /i = 2 in eq. (2.12) and color factor Tr^^l in eq. (2.19) 
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For /i = 3 in eq. (2.13) and color factor Tr^^l in eq. (|230|): 
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For h = A in eq. (2.14) and color factor Tr^^' in eq. (2.20) 
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For /i = 5 in eq. ( 2.15 ) and color factor Tr'^' in eq. ( 2.21 ): 
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/i = 6 in eq. (2.16) and color factor Tr'^' in eq. (2.21) 
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For /i = 6 in eq. ( 2.16 ) and color factor Tr^^' in eq. ( 2.21 ): 
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B. Auxiliary functions for two-loop scheme shifts 
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